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Novel material properties can be realized by designing waves' dispersion relations in arti- 
ficial crystals. The crystal's structural length scales may range from nano- (light) up to 
centimeters (sound waves). Because of their emergent properties these materials are called 
metamaterials. Different to photonics, where the dielectric constant dominantly determines 
the index of refraction, in a ferromagnet the spin- wave index of refraction can be dramatically 
changed already by the magnetization direction. This allows a different flexibility in realiz- 
ing dynamic wave guides or spin- wave switches. The present review will give an introduction 
into the novel functionalities of spin- wave devices, concepts for spin- wave based computing 
and magnonic crystals. The parameters of the magnetic metamaterials are adjusted to the 
spin- wave /c- vector such that the magnonic band structure is designed. However, already the 
elementary building block of an antidot lattice, the singular hole, owns a strongly varying 
internal potential determined by its magnetic dipole fleld and a localization of spin-wave 
modes. Photo-magnonics reveal a way to investigate the control over the interplay between 
localization and delocalization of the spin- wave modes using femtosecond lasers, which is 
a major focus of this review. We will discuss the crucial parameters to realize free Bloch 
states and how, by contrast, a controlled localization might allow to gradually turn on and 
manipulate spin-wave interactions in spin-wave based devices in the future. 
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30 1. Introduction 

Magnonics is a young and evolving research field. Its aim is to control and manipulate 
spin waves in ferromagnetic materials [Ij. In analogy to photonics taking control of the fiow 
of light, it allows for the design of material properties for spin waves. The experimental 
realization of computing with spin waves will be discussed first. These are logic gates 

35 (NOR, XNOR, etc.) - novel devices that benefit from the particular properties of spin 
waves. To connect and to guide information, spin-wave frequency filters and wave guides 
are under research. In such, novel computing concepts have been described. One of these 
are reconfigurable mesh structures using spin waves for parallel graphics computing. This 
review will introduce how periodic structures can be realized in diflFerent dimensions, which 

40 will determine the line, plane or volume in which the artificial material properties are tailored 
and affect propagating spin waves. A large part of the present review will map how they can 
be investigated using laser light ^ hence the term photo-magnonics - a technique developed 
in the last years. In particular, two-dimensional structures will be considered. Generally, the 
response function to a dynamic excitation propagating through a periodic crystal determines 

45 its propagation defined by the dispersion E{k)^ which in turn determines the available density 
of states. For weak potentials, the concept of constructing Bloch states using plane waves 
and calculating their dispersion in the reduced zone scheme is straight-forward. Applied 
to photons in periodic dielectrics, we show how band gaps at the Brillouin zone boundary 
tailor novel functionalities. 

50 In the case of spin waves the dispersion can be quite complex in the unstructured mate- 
rials already; we give a short introduction into the peculiarities of the spin-wave dispersion, 
namely the anisotropy of the dispersion with respect to the magnetization direction. The 
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same concepts as for electrons and photons can be applied to magnons to form Bloch states 
in a magnonic crystal. In the uniform mode approximation simple band structures can be 

55 directly calculated. However, these ignore the detailed shape of the potential and are a good 
approximation for a weakly varying potential only. A more realistic potential landscape will 
be discussed in the case of the periodic dipolar potential for a simple antidot lattice. The 
individual shape of each spin-wave scatterer, the building block of the magnonic crystal, 
is determined by its detailed variation of the dipole field. Spin-wave modes in periodically 

60 structured materials have been studied experimentally by different techniques: microwave- 
based techniques (ferromagnetic resonance (FMR), vector network analyzer FMR or pulse- 
inductive microwave magnetometer (PIMM)), in combination with spatial-resolved Kerr 
microscopy. A somewhat larger part will be devoted to optical pump-probe techniques with 
femtosecond lasers in this review. The interplay between localization and delocalization of 

65 the dynamic modes in the magnonic crystal is shown to be crucial to understand the observed 
spin- wave modes. The deep distortion of the internal field by the dipole field at an antidot 
site results in a localization of the spin waves in many cases, prohibiting the observation of 
delocalized modes. However, in order to realize magnonic wave guides and active spin-wave 
materials, it is crucial to control the interplay between localization and delocalization of the 

70 dynamic modes in the magnonic crystal. For a high filling fraction of the antidot lattice 
(i.e. amount of material removed), spin waves are trapped over a large field range at places 
with reduced internal field in the magnonic crystal. Generally, only for small filling fractions 
delocalized free Bloch modes are observed, whose /c-vector is determined by the crystal's 
unit cell. We will discuss how these Bloch modes can be identified in photo-magnonics. At 

75 the end we will give a short outlook on future possibilities, perspectives and developments 
in the field, such as guiding spin waves in magnonic wave guides, spin- wave resonators or 
tuning the degree of localization of a spin wave and their interaction. 

2. Computing with spin waves 

On the semiconductor roadmap of 2009, possible emerging research devices and emerging 

80 materials for future electronics and logic are presented [2j. On the materials side, graphene, 
carbon nanotubes, and nanowires are discussed as possible candidates to develop materials 
to meet requirements for future semiconductor devices. On the device side, the use of the 
electron spin for information storage, and the implementation of spin waves for performing 
non- volatile logic functions on a CMOS chip are discussed. The idea to use waves for signal 

85 processing is not new: light in glass fibers is used as an interconnect between chips in 
powerful computers. The IBM blue gene super computer link modules have optical printed 
circuits on each board. However, their function is to convert electrical signals into optical 
signals that can be directed to the next step of data processing. As long as the processing 
relies on electrons the extra effort of transforming the data into light pulses for effective 

90 and fast throughput to other regions of the processor is unavoidable. Building a computer 
that uses solely photons for information processing and its non-volatile logic functions is 
not that easy a task. It is difficult to realize a light switch that can be operated by light; 
to switch light, one would need to modify the index of refraction, which is related to the 
electron density of the material. Consequently, a high density of photons is needed to have 

95 a significant effect on the dielectric properties of a material. To realize a switch that should 
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turn on and off the total number of photons, certain tricks have to be played. We will discuss 



How is it possible to take advantage of the properties of spin waves, the intrinsic strong 
non-linearities to realize switches for spin waves and to realize non-charge based persistent 

100 devices for spin-wave based computing? In recent years a large portion of the effort of the 
magnetism community was devoted to the field of spin electronics, to develop spin-based 
semiconductor devices for spin-based transistors and spin-based logic. Magnetic random 
access memories (MRAMs) have been developed that exploit the advantage of non- volatile 
memory storage, while reducing power consumption. Power consumption is one of the 

105 biggest problems of high integration circuits today, because of the heating of the chip, re- 
duced battery lifetimes and environmental harmfulness. The idea of looking at spin waves 
stands to reason, which has, however, not been done until recently. Even though spin waves 
have been studied as narrow frequency filters in the 60s to 80s in YIG-based high frequency 
devices [3^ 4J, the idea to use spin waves for data transmission and processing is very new. It 

no can be realized with control on the nanometer scale by now and will be discussed in section [5} 
What makes spin waves favorable for technological applications is the corresponding index 
of refraction, which can easily be manipulated. The term 'index of refraction' should here be 
understood in analogy to photonic crystals where optical resonances can be shifted by a pe- 
riodic modulation of the dielectric constant. In the ferromagnetic case, the modulation may 

115 be in the parameters exchange constant (^4) and/or saturation magnetization (Ms) while an 
additional directional asymmetry is introduced by the applied magnetic field. Corresponding 
Bloch-like resonances of the spin- wave spectrum will be discussed in section [KS) In general, 
spin-wave frequency and wave length can both be tuned from MHz to THz on micron to 
sub-nm scales, respectively. Therefore, the use of spin waves for logic devices is the natural 

120 extension of magnetic non-volatile elements for storage. The flexibility of spin-wave based 
materials opens up the possibility of a spin- wave based computing architecture. This is the 
underlying idea of magnonics: to use spin waves for power-saving computing. Nevertheless, 
certain prerequisites in manipulating and guiding spin waves have to be demonstrated a pri- 
ori [5^£j: (i) the controlled excitation of spin waves at defined frequencies, (ii) guiding spin 

125 waves in magnetic wave guides and (iii) active spin- wave devices for spin- wave manipulation 
and information processing. 

2.1. What are spin-wave guides? 

An object often referred to as a spin- wave guide is a (quasi) one-dimensional stripe of 
ferromagnetic material. Through its physical boundaries, it restricts the region of existence 

130 for spin waves. In optics, an analogue object is called a fiber. Different from glass fibers in 
optics, a spin- wave guide is a homogeneously magnetized wire that can channelize, split, and 
manipulate submicrometer- width spin- wave beams, as was shown by Demidov and Demokri- 
tov [5j. The differences are depicted schematically in Fig. [l] where the properties of a graded 
index glass fiber are compared to a magnetic wave guide. In the graded index glass fiber, the 

135 index of refraction gradually increases to the outer shell of the glass fiber having a higher 
index. The light is totally reflected and conducted in the core of the fiber. The gradual 
transition leads to a smooth curvature of the wave fronts so that the Gaussian shape of 
the intensity profile and temporal coherence of the signal pulse is only slightly distorted. 
The index of refraction, arising from the dielectric properties of the materials, serves as a 
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Figure 1: A graded index glass fiber (left column) concentrates the intensity of the light modes in the 
area of low index of refraction, i.e. small dielectric constant (depicted by a solid black line). In a simple, 
homogeneously magnetized plane magnetic wire (center column), the internal effective field (solid black 
line) determines the localization of the spin- wave modes. Localization of spin waves along a spin- wave guide 
(right column) measured with micro-BLS (Brillouin Light Scattering) (adapted from ^). The internal field 
changes as the width of the wire changes, resulting in different localization axes of the beam path. 



potential landscape for the light wave. For a spin wave, the situation is sketched in the 
middle of Fig. [T} A flat, structured stripe with a width in the micron range serves as a spin- 
wave guide. The details of propagation depend on the magnetization direction: magnetized 
in-plane perpendicular to the wire axis, Damon-Eshbach modes that are a species of dipolar 



spin waves (see section 4.1) can travel along the wire. Magnetic charges at the boundaries 
lead to a decrease of the internal magnetic field in the wire at both sides; i.e., due to the 
demagnetization field of the wire magnetized perpendicular to the wire axis, the effective 
field is reduced at the edges of the spin- wave guide (see sketch in Fig. [T]). The well structure 
of the internal field shows minima at both sides of the wave guide. Depending on the width 
of the stripe, the internal magnetic field in the center of the stripe is also diminished more or 
less in strength for smaller or wider wave guides. One can excite modes, which are localized 
in these wells ^ so-called 'edge modes', as well as modes in the middle of the wave guide - 
so-called 'center modes'. In this context, it is important to notice that edge-modes always 
have lower frequencies, as compared to center modes. 

The right panel in Fig. [l] shows micro-Brillouin Light Scattering (micro-BLS) maps of 
such modes. Spin waves are excited in the region with smaller width, and thus stronger 
reduced internal field at the center. They propagate towards the region with increased 
width and thus increased internal field in the center. A transformation from a center into 
an edge-mode takes places in the transition zone. Due to the upward shift of the spectrum 
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in the wider stripe, the excited wave can no longer exist in the center. This locahzation 
160 process depends on the frequency of the initiaUy excited mode: for higher frequencies, the 
edge-modes move toward the center. This example nicely shows that the understanding of 
and control over the spin-wave dispersion, localization, and delocalization with the internal 
field distribution is one of the key aspects to progress in this novel sub-field. 

2.2. Spin waves on a chip: Reconfigurable mesh design 

165 Concepts to use spin waves for data processing in a chip have been developed by Kang 
Wang et al., a short review on spin- wave based computing has been published in the series 
of the first "International Seminar and Workshop Magnonics: From Fundamentals to Appli- 
cations" held August 2009 in Dresden [6] based on their earlier work published in [7] . Many 
current applications require a vast amount of data to be processed in parallel. A typical 

170 example is a graphic chip, with its algorithms for parallel image processing. Their architec- 
ture and computation power has been increased enormously in the last decade. From the 
scientific point of view, numerical simulation programs, which can be parallelized, can take 
advantage of the significant progress in that field, so that when put together, an array of 
graphic chips makes a new super computer. The idea of the spin- wave based computing con- 

175 cept to parallelize computation uses a reconfigurable mesh architecture. Spin- wave guides 
transmit the signal at each line. The chip consists of a mesh of x spin- wave switches 
interconnected by ferromagnetic spin- wave guides. Each node is realized by a ferromagnetic 
switch. If the switch is "on" , the spin wave is guided into the crossed line to the spin- wave 
buses output. The switching frequency is in the order of GHz and transmission speed is 

180 lO^ms"^, allowing fast data processing. In the architecture presented in Refs. P,[2j, excita- 
tion is realized by a strip-line; detection will be realized by inductive detection at a second 
strip-line. The spin-wave switch at each crossing could be realized by a diluted magnetic 
semiconductor which can be switched from a ferromagnetic to a paramagnetic state by ap- 
plying positive or negative voltages. It should be remarked that as of today, this concept is 

185 lacking a realistic practical implementation. Ferromagnetic semiconductors like Mn-doped 
GaAs have both a very strong spin-orbit interaction and a critical temperature that is be- 
low room temperature, and hence will not work for actual devices. However, using a mesh 
structure allows for large image data to be processed in the reconfigurable nodes. This is a 
typical application needed for example to process image data for compression in a television 

190 graphics chip. Additionally, different to a standard reconfigurable mesh design, spin-wave 
buses could in principle address different frequencies. Therefore, parallel operation seems 
possible as well. The layout is shown in Fig. [2j In a), the schematic layout of one node is 
shown, which can be addressed by the voltage turning on/off the magnetization, which then 
guides the spin- wave package to the lower magnetic wave guide. In b), the full layout of the 

195 mesh structure is depicted. 

2.3. Elements of spin-wave computing and logic 

At the same time as concepts that take advantage of spin-wave computing have been put 
forward, spin- wave based logic elements have been demonstrated. Their principle relies on 
a spin-wave beam splitter. Two spin-wave packets of same amplitude are split into different 
200 arms of a spin-wave guide. In one arm, the propagating spin wave receives a phase shift 
by 180° (tt). The arms are united and both spin waves interfere destructively. This can 
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Figure 2: Concepts and first realization using spin- wave based computation: a) Schematics of a spin- wave 
bus, and b) building blocks of a reconfigurable mesh structure (from [7 ). c) Experimental realization of 
a XNOR functionality (logic equality) using a spin-wave interferometer with phase shifter as reproduced 
from |9j. 
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be compared to a spin wave Mach-Zehnder-type interferometer. Different realizations for 
phase shifters have been suggested; for example, in micromagnetic simulations, it was shown 
that a domain wall could be utilized as a phase shifter in one of the interferometer arms [8j . 

205 More generally, the phase shift, slowing down or acceleration in one interferometer arm 
relative to the other, depends on the dispersion of the spin wave. The dispersion is naturally 
modified by the (sometimes complicated) magnetic structure of the domain walls, thus the 
average internal field in the domain wall. Another possibility is to shift the frequency by 
locally applying an external magnetic field. In the yttrium iron garnet (YIG) wave guides, 

210 which have the lowest damping of any magnetic material, millimeter propagation length of 
spin waves in the microwave frequency range can be realized. In these systems, spin- wave 
propagation, spin-wave amplification, spin-wave pumping, phase shifters, interferometers 
and filters were demonstrated [91 Schneider and coworkers developed in these wave 

guides a spin-wave logic realizing exclusive-not-OR and not- AND gates based on the Mach- 

215 Zehnder-type interferometer, which will be presented in more detail. Input and output into 
the YIG wave guides is implemented by microwave antennas; the signal is propagated in 
the YIG film and split into the two arms of the interferometer. The phase shifter is realized 
by applying an Oersted field on top of the YIG, which changes the carrier wave number of 
the spin- wave packet within that region. By changing the amplitude of the Oersted field, 

220 the phase shift can be chosen to sum up to one half of a wavelength. The output signal is 
shown in Fig. [2]c). If no field is applied at either arm, both spin- wave packets are equal 
in phase and the full signal is detected. If for any of the arms the relative phase is shifted 
(input 0,1; 1,0), the signals detected at the microwave detection antenna have opposite signs 
and cancel out. The output signal is zero. If both arms are subject to an Oersted field, 

225 both spin waves are shifted by the same phase, and again, the full signal is measured at the 
output. This gives the functionality of an XNOR logic device (logic equality). A second 
logic device functionality is implemented by using the Oersted field on top of the YIG arm 
to completely suppress the propagation (spin- wave switch) , which gives a zero transmission 
(output 0) if both inputs are on (input 1,1), else (input 0,0; 0,1; 1,0) a signal is detected 

230 (output 1), which is the NAND functionality. The disadvantage of YIG-based magnonic 
devices is that YIG cannot be integrated into standard semiconductor technology because 
high quality YIG films cannot be grown on silicon at the moment, but special substrates 
are needed. The industrial demand of miniaturization and integration into semiconductor 
technology is much better met by Permalloy (Ni8oFe2o). This is why the aforementioned 

235 concepts need to be transferred to smaller structures. 

In Fig. |3] Wang et al. compare MOS logic and magnonic logic in terms of throughput (the 
number of operations per area per time) as a function of the minimum feature size A, which 
is the gate length for CMOS or the wavelength for a spin- wave circuit, respectively. The 
throughput in devices realized so far is quite small and feature sizes lie in the micrometer 

240 range. However, according to the projected estimates, spin logic may provide a throughput 
advantage of more than three orders of magnitude over CMOS in the future. Magnonic logics 
have several scaling advantages, one being the inversely proportional throughput of the spin 
circuit to the wavelength. However, scaling down to spin- wave length in the nanometer 
range is in first order a material problem: one has to find materials where the damping 

245 of spin waves is low [llj. Half metals are an ideal candidate here, since electron and spin 
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Figure 3: CMOS technology scaling versus spin- wave based magnetic circuits, from [6 . For spin- wave 
circuits, the number of operations per area per time (throughput) scales inversely proportional with the 
minimum feature size (wave length) A. 

scattering channels can be decoupled [12]. 

3. Tailoring artificial materials 

To realize devices described in the previous section, it is most convenient to microscop- 
ically manipulate the properties of a single material. This is the origin of the success of 

250 the semiconducting materials as silicon, which has a band gap at the Fermi level. Thus, 
no states are available for conduction inside the gap region. To design the electronic prop- 
erties, doping is necessary. The potential landscape is then formed by a combination of 
differently doped regions and various devices such as transistors can be realized to switch 
current flow through a device on and oflF. It is the flexibility of the material that allows a 

255 change in the conductivity by orders of magnitude and allows the shaping of the density of 
states at well-deflned energies with a spatial precision down to nanometers. The band gap 
and a controlled introduction of defect states deflnes this flexibility. In a similar manner, 
photonic materials can be designed. Spectacular examples are materials with negative in- 
dex of refraction [T3j, photonic band gap materials already existing in nature giving it its 

260 colorful life, and controlled "defects" realizing high quality resonators for light waves. These 
materials own novel properties that are very diflFerent from those of their constituents. Any 
wave propagating probes a material with effective parameters designed by structure and 
composition. Therefore metamaterials are per deflnition artiflcial materials where the pe- 
riodicity of the structure is smaller or equal than the wavelength. This can be centimeter 

265 or millimeter in size: flrst metamaterials have been realized for microwave radiation [14j. 
Artiflcial materials also include those acting on water waves [15j, which allow the formation 
of refractive elements like lenses and sonic waves, to realize sonic crystals with negative re- 
fraction [16j. These examples from other areas utilize rather macroscopic waves transmitted 
through the metamaterial. In the following, we will remind the reader of the basics of band 

270 structure formation in metals and covalent bonded materials. As a next step we will apply 
this formalism to photons in a periodically structured material with alternating index of 
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Figure 4: Schematic analogy of electrons, photons, and magnons in periodic potentials. The dispersion 
is shown in the reduced zone scheme. Different dispersions result in characteristic differences in the band 
structure. A common feature is the possibility of the formation of band gaps. The polarization (transverse 
electric (TE) or transverse magnetic (TM)) for photons and the propagation direction /c, with respect to the 
applied magnetic field for magnons (parallel or perpendicular), are additional free parameters determining 
the band formation for the latter, (/c-scale is logarithmic for magnons) 

refraction before we expand the formalism to magnons. 

3.1. Block conditions and hand structure: weak periodic potential 

In a metal, the free electron ansatz is very successful, leading to Sommerfeld equations 
describing transport in an electron gas. The kinetic energy dominates the properties of the 
electrons. The energy levels are filled homogeneously up to the Fermi level leading to a 
Fermi sphere in momentum space, which behaves isotropic for the different directions. The 
Hamilton operator is given by 

The solution for the free electrons where V{r) is neglected, is simply given by a quadratic 
dispersion mirroring the increase of kinetic energy 

2m 
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Figure 5: Free electron band structure (left) in the periodic boundary scheme and reduced zone scheme. 
In a weak periodic potential a band gap opens up (middle). It is given by 21^01 where Vg is the Fourier 
component of the periodic potential. Free electron band structure (right) for a fcc-lattice along the main 
symmetry directions (adapted from [17j). 

Turning on a periodic potential, the Bloch theorem allows to exploit the periodicity. The so- 
lution can be separated into two parts: a free electron propagating part with wave number k 
and a function that is periodic with the lattice Un^k- The new quantum state is described 
by a wave number k and a band index n, 

V = V{r + R), 

The periodic function Un^k describes the electron distribution within one cell around the 

275 crystal atom. Without making any calculations of the exact band structure, already some 
general conclusions regarding the electron states in the resulting band structure can be 
drawn. As a consequence of the periodic potential, the electron dispersion mirrors the 
symmetry of the lattice. Gaps open up especially around the zone boundary in the high 
symmetry directions. This can be seen in Fig. [5] (middle) where, for a weak potential, 

280 the opening of the gap at the zone boundary is given. The periodic functions Un^k and 
their eigen energy En^k have to be calculated. To do this for a given potential, the plane 
wave ansatz exp(i/cr) is used to make a Fourier expansion for the wave function. It is 
most convenient to solve the set of Bloch states in fc-space via the Fourier transform of the 
Schrodinger equation. Then the Fourier coefficients of the periodic potential V determine 

285 the coefficients of the Fourier expansion of the wave function. This directly allows the 
calculation of the band gap: it is given by 2\Vg\ where Vg is the Fourier component of the 
periodic potential at the zone boundary k = G. In fact, we will see in Chapter 4 that for 
spin-wave states, similar approaches to calculate their band structure can be used, however, 
also with distinct differences arising form the different equations of motion including the 

290 symmetry of dipolar interactions mirrored by the vector nature of the magnetization M. 

3.2. Electrons in a strong periodic potential 

If the electron wave function possesses a localized nature, the free electron approach is 
not successful. The electronic band structure of such a material mirrors more strongly the 
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orbital nature of the atomic wave functions. The s,p,d-hke symmetry determines many of 

295 the fundamental properties of the material. Examples are band ferromagnets, where the flat 
d-like electron bands are responsible for the formation of the ferromagnetic order [18] . Other 
examples include carbon, silicon and germanium of the group 4 elements in the diamond 
structure. The local electron density is described by the linear combination of s and p wave 
functions forming the sp^ tetraeder structure. In that case, another possibility is to take 

300 the atomic wave function of the single atom as a basis to construct the bands. Typically, a 
Wannier function localized on an atom, but with decaying density of states to the neighboring 
atom, is used. The band width is given by the overlap to the next atom. A large overlap 
means a large hopping rate and a strong delocalization and consequently a large bandwidth. 
The more localized the electrons are on the atomic-like orbitals, the smaller is the overlap 

305 and the smaller is the band width in energy. One can derive these tight binding bands 
for different crystal symmetries (Fig. |6]for bcc and fee). For example, in the bcc case, the 
splitting into the dxy, d^zj dy^ symmetry and dsz'2_j.2^ d^2_y2 symmetry are observed varying 
from a different distance to the next neighbor atom along the axes of the cube and along the 
diagonal. The band width of the d^z2_j.2^ d^2_y2 bands is increased by a factor of 1.5 along 

310 the F-X symmetry direction. The symmetry of the atomic wave function in connection with 
the lattice symmetry determines the electron states observed in the crystal; even though the 
states are strongly determined by their atomic wave functions and show a weaker dispersion, 
the next neighbor distance significantly infiuences the properties of the electron in the band 
such as localization. We will see here the similarities to the magnonic systems, where spin 

315 waves can also be of localized nature later and its symmetry is mirrored by a short range 
interaction. 
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Figure 6: Canonical tight binding d-band structure for fee and bee lattiee along the main symmetry diree- 
tions. The band width ujd of the eosine-like bands is given by the overlap td of the atomic orbital functions 
(here d-orbitals) into the different crystal directions ~ 2^Nnntd' The overlap td determines the delocaliza- 
tion of the electrons within the band. The symmetry is determined by the underling atomic wave functions 
admixture (adapted from [171 HE] ) • 
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3.3. Photonic crystals - photons in periodic potentials 

Similar to the periodic potential experienced by electrons around ion cores in a lattice, in 
a photonic crystal the periodic dielectric properties determine the periodic confinement for 
the photons. In general, if the potential varies periodically for a wave traveling through a 
two-dimensional artificial crystal, one can obtain a solution by labeling the modes with their 
band index n and wave vector k using Bloch's theorem, as described for the electrons before. 
The same concepts can be applied despite the equation of motion being determined by the 
wave equation derived from the Maxwell equations. Equal to the periodic potential of the 
positive ion cores Vu in the solid state for electrons, for a photonic material the dielectric 
properties e determine the refractive index, which determines the periodic confinement for 
the photons [19j. To compare: in a magnonic material, the dispersion will be dependent on 
the local magnetization direction Vm and the induced internal field. In many cases photonic 
structures are two dimensional. The mathematical ansatz to describe the electro-magnetic 
modes of a photonic crystal is similar to the description of delocalized electronic states in 
solids. The wave traveling through a two-dimensional artificial photonic crystal experiences 
a periodically varying potential arising from the dielectric constant e. In two dimensions 
the potential varies as 

Ve = V,{r + R), 

where i? is a linear combination of the primitive lattice vectors for a two dimensional array 
in the x-y plane. Using Bloch's theorem and using p as the projection of r in the x-y plane, 
one obtains a solution labeling the modes with their band index n and wave vector k\\ 

One example taken from [19j is given in Fig. [7| It has been calculated for a two-dimensional 
array of aluminum columns with the dielectric constants of Al [e = 8) and air = 1). 

320 In the case of an electro-magnetic wave, one has to distinguish two solutions in the two 
dimensional crystal. For a photon propagating through the crystal, the electric field can 
be polarized in the plane or out of the plane of the two dimensional structure. Therefore, 
different dispersions are found for transverse electric (TE) and transverse magnetic (TM) 
polarization to the plane of the crystal: TE and TM modes show a different slope at the F- 

325 point. The TM-mode shows a lower velocity because it effectively sees the higher dielectric 
constant of the Al columns. The corresponding modes at the zone boundary along the 
different symmetry direction at the X and M point are plotted to the right: they form 
standing waves which are localized either at the Al columns or in between. The one with 
lowered energy is localized in the areas with high s and the one with increased energy in the 

330 areas with low 5, leading to distinct gaps in the photonic band structure in this direction. 
It can be easily understood that the gap is larger for the TM modes for the same argument 
that this mode, since its £'-vector is in the plane, is strongly affected by the alternating 
dielectric constants. 

3.4' Novel functionalities of photonic crystals: slow photons 

335 A designed photonic band structure can be used to achieve novel optical properties. 
Quite generally, close to the zone boundary, the band fiattens out. Approaching the zone 
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Figure 7: Photonic band structure calculated for an array of two dimensional Al columns. Transversal 
electric (TE) and transverse magnetic (TM) modes show a different slope at the F-point. The TM-mode 
shows a lower velocity because it effectively sees the higher dielectric constant of the Al columns. Shaded 
areas mark the splitted bands at the high- symmetry points point. Right: At the zone boundary for the 
different bands, the standing waves with high and low energy are plotted (Al {e = 8.9) in air {e = 1)) 
(adapted from |19] ) . 

boundary then means in terms of group velocity given by duj/dk^ that the photon is "slowed 
down" [20j. Usually switching light by light is not possible: the change in the index of 
refraction is only small. However, when using slow hght, a trick can be applied. Exploiting 

340 the sensible reaction to the refractive index in the region where the band flattens out, it was 
experimentally demonstrated that it's feasible to guide light from one wave guide to another 
output, which is the realization of an ultrafast nanophotonic switch using slow light [21j. 
In detail, the spatial beating due to odd and even super-modes in two parallel wave guides 
is used for rerouting triggered by a femtosecond laser pulse. This shows that a wave guide 

345 in a photonic material is not only a straight line of material conducting light as a glass 
fiber: the peculiar properties can be designed by the photonic structure around the wave 
guide giving the light completely new properties in the artificially structured material. We 
suggest it should be one aim in magnonics to put forward similar concepts to exploit the 
full possibilities of artificially structured materials for spin waves. 

350 4. Magnonic crystals — spin waves in periodic potentials 

We will now in more detail discuss the effect of structuring a magnetic film to form a 
magnonic crystal. That is, only two-dimensional systems are taken into account and be- 
ginning with a continuous film, the effects of a periodic structures are considered on an 
analytical as well as numerical level. We will show that for the dipolar modes, the disper- 
355 sion differs for different propagation directions k relative to the magnetization vector M. 
Namely, the Damon-Eshbach {k _L M) and the backward volume mode [k \\ M) exist, which 
are different due to dipolar interactions. As a consequence, the spin-wave states in a square 
magnonic crystal have a lower symmetry than the underlying lattice structure. The manip- 
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ulation of the local magnetization direction or an applied field on a nanometer scale allows 
active dynamic control of the spin-wave diffraction on the nanometer scale. Thus, a rotation 
of the magnetization relative to the lattice changes the energy landscape dramatically. 

^.7. Spin-wave dispersion from nanometer to micron range 

Before we discuss the propagation of spin waves in a structured medium, in analogy to the 
preceding chapters on electrons and photons in periodic potentials, we first have to say some 
words to the equation of motion for spin waves to be solved in the periodic environment. In 
a continuous, non-periodic system, the magnetization follows the Landau-Lifshitz-Gilbert 
equation of motion, which reads 

-7/xoM X i/eff + — M X . (1) 



dt Ms V dt 

We will only consider a field applied in the plane of the film. The case of perpendicular 
magnetization is technologically less favorable since it demands strong fields for saturation, 
and it is also physically less interesting because the spectrum is isotropic. In the macrospin 
approximation, all the individual spins are considered to precess in phase and the thin-film 
solution to Eq. ([T]) is given by the Kittel equation 

'''^' = hJh^ + Ms-^). (2) 



7Mo/ V l^oMs 

Here, represents an effective anisotropy in the out-of-plane direction. If one allows for 
solutions to Eq. ([T]) other than uniform precession, it becomes important to distinguish 
two different interactions coupling magnetic moments or individual spins, respectively. The 
corresponding solutions will be discussed in the following sections |4.1.1| and |4.1.2[ 



4.1.1. Dipolar spin waves 

Solving Eq. ([T]) in the magnetostatic limit yields a manifold of dynamic solutions for 
a thin magnetic film of thickness t [23l EH [25j . All can be classified by their propagation 
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Figure 9: Schematic spin- wave dispersion. In the micron wave length (left), dipolar interactions dominate. 
The Damon-Eshbach and backward volume modes are shown. At an arbitrary angle in between further 
modes are found (gray shaded area). In the nanometer region (right), the dispersion shows cosine-like 
behavior. Here, the exchange interaction is dominant and can be approximated by a parabola for small 
energies. It intersects in a broad region of high-energy spin- waves excitation, where spin waves are heavily 
damped (red shaded area). 

380 direction with respect to the magnetization. Modes whose frequency hes above the one of 
the uniform precession generaUy tend to locahze at the surface of the film and have a wave 
vector pointing perpendicular to the magnetization. The geometry with A; _L M || H is 
identified as the Damon-Eshbach (DE) geometry [24j. The particular relation between the 
magnetic field and mode frequency is given by 



385 For directions k \\ the so-called backward volume waves occur, with a reduced preces- 
sion frequency compared to the uniform precession. This leads to a negative dispersion, 
as the waves travel "backward" in phase. In between, one finds a manifold of spin waves 
corresponding to the continuous change of angle from parallel to perpendicular (grey shaded 
region in Fig. [9]). This peculiar energy dependence on the angle arises from the dipolar 

390 interactions. For wave lengths below a micron, the energy difference becomes smaller and 
both dispersions are degenerate in energy. Here, the exchange interaction becomes impor- 
tant, so that this contribution has to be taken into account for mixed dipole and exchange 
spin waves in an intermediate region of length scales [23j . The steep uprise of the energy 
for smaller wave length marks the dominance of the exchange interaction. In a thin film, 

395 the confinement perpendicular to its plane is dominated by the exchange interaction, while 
in lateral directions no such restrictions exist and dipolar magnetostatic spin waves may be 
formed. 

4.1.2. Exchange spin waves 

Since spin-wave lengths span several orders of magnitude from tens of microns (even 
400 higher for low-damping materials) to below 1 nm, also their frequencies may vary from GHz 




(3) 



to THz. In addition, the frequency for a given wave length can be shifted by the magnetic 
field. This broad region in length and time scales is one reason that makes spin waves so 
interesting for high frequency applications. However, also the dominating interaction varies: 
at wave lengths below 100 nm, the dispersion is dominated by the exchange interaction. 

405 The magnetostatic contribution to the energy of the wave can be neglected. This simple 
picture is solved in many solid states physics textbooks for a chain of precessing spins, where 
next neighbors are coupled by the exchange interaction. The solution is a (1 — cos(/cr)) like 
behavior. As a consequence of neglecting the anisotropic dipolar contribution, the dispersion 
in the exchange limit does not change with the magnetization direction. It only depends 

410 on the next neighbor distance r and the strength of the exchange interaction, and can be 
calculated in the "frozen magnon" picture from the electronic structure. For small /c, the 
dispersion can be approximated to a quadratic form. In that region, the energy increases 
quadratically with momentum hk similar to the free electron behavior. One can think of 
realizing "free electron like" magnonic materials. 
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Figure 10: Schematic spin- wave dispersion for the different length scales in the periodic zone scheme. For 
a spin- wave Bloch state, prominent effects in the band structure are expected: 10 /im range (left), 1 jam 
(middle) and 100 nm (right). 

At even higher energies, the spin waves approach the THz regime. Their energy gets 
comparable to single spin-flip excitations (Stoner excitations) between the bands of diflFerent 
spin character. These high-energy spin-wave modes in the THz range are heavily damped, 
indicated by the red shaded area. Dynamic excitations have to be thought of as a superpo- 
sition of multiple spin excitations propagating through the ferromagnet [261 EZ] • Life times 
of these high-energy spin waves of nm length scales are a few picoseconds before the decay 
into other spin- wave excitations of lower energy [28j. Investigated heavily in the 80s |29j . 
their physics have drawn new interest. They are thought to be the key to the further under- 
standing of ultrafast demagnetization processes after femtosecond laser excitations. Because 
of their short wave length and lifetime, at the moment, these are not suited for studies of 
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425 the formation of magnonic bands. They propagate only a few nanometers. In the foUowing, 
we wiU discuss the formation of magnonic band structures in the dipolar region. 



4^2. Band structures of magnonic crystals 

In this section we will discuss first what are the expectations for a band structure in 
a magnetic material that is periodically modified. If novel spin- wave states appear due to 
the periodic modification, a magnonic crystal is formed. Second, we will calculate for one 
example the respective band structure. For wave lengths larger than 1 /xm, the dispersion of 
the spin waves is dominated by dipolar interactions. Corresponding frequencies are below 
20 GHz, depending on the magnetic material and applied field. A steep quadratic increase is 
related to the nanometer wave length range. The strong magnetic exchange interaction takes 
over as the frequency reaches 1 THz. The dispersion has completely different slopes for both 



regions. The consequence is shown in Fig. 10: in the short wave length range a periodic 



modification will lead to a "quasi free electron" -like band structure in the periodic zone 
scheme (schematically constructed on the left), similar to what is found in micromagnetic 
simulations [30]. If dipolar interactions dominate, the energy splits for the Damon-Eshbach 
and backward volume modes: the band structure will be anisotropic with regard to the 
applied field. When the structure size is larger than the micrometer range, the negative 
dispersion of the backward volume mode leads to the surprising result that the second band 
can be lower in energy than the first band. We will see that this unexpected feature is found 
as well in our band structure calculations. A very complex spin- wave manifold can develop 
with many bands having similar energies. The three types of magnonic crystals classified in 



Fig. [To] will show a very diflFerent behavior: dipolar bands with negative or positive dispersion 
as well as quadratic exchange dominated bands. 

In the following, we will discuss in greater detail how, similar to the Bloch theorem 
applied to electrons and photons, a magnonic band structure of a periodic ferromagnetic 
system can be computed. A theory was developed by Puszkarski et al. ^TJ [S^- A periodic 
alteration of the saturation magnetization was put into the Landau-Lifshitz Gilbert equa- 
tion and solved by a plane-wave method. This theory is based on infinitely extended, three 
dimensional magnonic crystals. While this ansatz emphasizes the close relationship of pho- 
tonic and magnonic crystals, it neglects the particular property of inhomogeneous internal 
magnetic fields by setting Hi^t = i/gxt- Note, that only if the regions with an inhomoge- 
neous i/int are restricted to a vanishing fraction of the magnonic crystal's unit cell, this can 
be justified. This is for example the case for a two-dimensional magnonic crystals, whose 
typical structure size is much smaller than its thickness. We will in the following show, how 
the plane-wave ansatz can be modified to apply to such a magnonic medium - exemplarily 
represented by an antidot lattice. 

Spin waves with a wave length A which is much bigger than the film's thickness t (external 
field in the film's plane), have almost uniform mode profiles in the z-direction. Thus, in the 
following the lower order modes' profiles will be assumed to be uniform and the so-called 
uniform-mode-analysis will be applied (see Hurben and Patton for details [33j). The basic 
idea of this approach is to neglect the exchange interaction, valid for sufficiently small k as 
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opposed by the full theory in Ref . [23] , and solving the linearized Landau-Lifshitz equation 



iQrriy — + -^hz = 
H 

Ms 

iVtrriz + rriy Jf^v ^ ^* 

The ansatz for the dynamic magnetization in the framework of the uniform-mode-analysis 
reads: 

rriz {x,y) = m^e'(^"^+^^^^ . (5) 

In the above equations rUy and are constant across the thickness. From the electromag- 
netic boundary conditions, it follows that the dynamic magnetic field reads: 

kt 

hy -my—e^^^''^^^yy^ sin^ 

By combining Eq. ^ and Eq. ^ into Eq. Q , the following direction-dependent dispersion 
relation can be derived: 

In Eq. Q H^y = H + Msf sin^, H^^ = i7 + Ms - Msf , where Ms is the saturation 
magnetization, g is the gyromagnetic ratio and (j) is the angle between the external field and 
the wave vector. As a next step, a periodic modulation of the saturation magnetization is 
introduced: 

Ms(r) = 5]Ms(G)e^^^ (8) 

G 

Where G = [G^, Gm]^ = ^^]^ is a two-dimensional vector of the reciprocal lattice. 

Any geometry can now be specified by an analytic expression for the Fourier components 
Ms(G). In order to circumvent additional boundary conditions arising from the air regions 
in the experimentally interesting geometry of antidot lattices, the following trick can be done: 
by filling the antidots with some artificial ferromagnet with a very high magnetic moment, 
spin waves in such a system either exist in the antidot, or in the surrounding matrix. Thus, 
a periodic energy landscape is constructed, which confines the spin waves. After solving the 
eigenvalue problem, one can remove the unphysical solutions, which predominantly dwell 
in the antidot. The rather smooth nature of the boundary, as a result of the cutoflF by the 
Fourier expansion ([s]) tends to mix the solutions for the matrix and the antidot. This results 
in non-vanishing imaginary parts for the frequencies of the modes in the matrix. However, 
one would need to use an infinite number of Fourier components to model a sharp transition 
between matrix and antidot, which would immediately violate the condition t. 
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In the case of magnetic discs (M|) with radius r, which are periodicaUy arranged on 
a square lattice with side length a = a' and embedded in a magnetic matrix (Mg), an 
analytical expression for the Fourier components of the magnetization profile can be found 
in ^ 

M.(G)^?^<M)/.(P), 

where P = r|G|, / = vrr^a"^ is the filling fraction and Ji is the Bessel function of first 
475 order. Note that a modulation of other material properties like the gyromagnetic ratio or 
the exchange stiflFness are not considered here. The generalization of Eq. ^ for the periodic 
material is a Bloch wave expansion: 

m{r) = ^mu[G)e'^^^^> (9) 

G 

Using the uniform-mode- analysis, the dynamic magnetic field components are 

G 

= Y: (-m,,.(G) '^+f^|-% inVc.) e^(^+^)^ 

G ^ ^ 



^^(k+G)r^ (10) 



where (j)G is the angle between the external field and K = k + G. The equations ([9]), (10) 
and Q provide a system of equations (finite number of lattice vectors) 

Mm{ = i^jrh{, (11) 

the eigen values of which have to be determined. Here, m-^ = [m^ ^(Gi), ^(Gat), ^.(Gi), m;^ ^[.(G 
The eigenvalues Qi = 271 fi {'yfiQH)~^ are proportional to the eigenfrequencies fi and the 
mode profiles can be constructed from the eigenvectors m-^. Note that the 2N x 2A^-matrix 
M has a block- diagonal form: 

In the equation above M^-^ = Mil = 
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with j = 1 . . . N . In this formula a locally varying static demagnetizing field i/dem was 
taken into account by including a Bloch wave formulation of this field, which is similar to 
the expression used for the magnetization profile: 



^dem(r) = V//dem(G)e^^^ (12) 



G 



In practice, the demagnetizing field can be obtained from numerical simulations, or from 
analytic expressions. Such a field may strongly alter the boundary conditions for the con- 
finement. Since the model for the antidots confines the waves in yet another way, it cannot 
be included for this approach. Nevertheless, artificial structures exist [34j, where the full 
model, taking into account the demagnetizing field and the exchange field, could similarily 
be applied. For a realistic situation, one problem has to be solved: the accurate description 
of a potential demands many reciprocal lattice vectors. But in this case the initial assump- 
tion X ^ t may be violated. A good agreement between theory and experiment may then 
only be expected for the lowest order branches. For the antidot model, we will only use 
thirteen lattice vectors here to discuss general trends. Note the sketched theory could be 
very easily altered to describe the situation, where - as for active devices - the spatially 
varying quantity is the applied field. 

Band structure calculations are performed for a CoFeB film with a thickness t = 50 nm. 



a hole distance a = 3.5 /xm and hole diameter d = 1 /xm (Fig. 11), and a nickel film with the 



same geometrical parameters (Fig. 13). A field of 130 mT is applied in the plane along the 
(1,0) direction. The discs are filled with an artificial ferromagnet, with a high saturation 
magnetization of 10 T. The artificial solutions, which predominantly dwell in the discs, are 
by virtue of the high magnetic moment shifted above the depicted frequency range. This 
ansatz separates the solutions into two sets. Only the low frequency which are of interest 



are given in Fig. 11 



Because of the larger saturation magnetization Ms, the bands for CoFeB are much higher 
in frequency. In the following we will describe some typical features. For comparison, the free 
spin-wave band structure is shown, with the first and second band marked. For this purpose, 
we have taken Eq. ([7| and plotted it in a reduced zone scheme. As expected, no band gaps 



at the zone boundary are present. On the left side of Fig. 11 the propagation direction is 



parallel to the applied field and thus determined by the backward volume geometry. As 



expected, with respect to the schematic bands given in Fig. [10] the second band is below the 
first band in energy. At the right, for the F ^ X' direction, the band dispersion is determined 
by the Damon-Eshbach geometry, showing a steep increase. The band calculation clearly 
reveals a different dispersion for the magnetization along the (l,0)-direction. Only for the 
M ^ F point the bands are the same. At the points of high symmetry, a splitting is 
observed which is of about 0.5 GHz in frequency. The uniqueness of the solutions inhibits 
sections of individual bands, as they appear in the free spin- wave picture. Instead, in the 
hybridization regions, the modes repel each other and interchange their character. Detailed 
theoretical studies of the group around Puszkarski showed the appearance of gaps in the 
band structure of two-dimensional and three-dimensional [32j materials where the spin- 
wave propagation is forbidden. They find that the size of the magnonic gap increases if the 
contrast between the constituents in magnetization or exchange stiffness is increased. 
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Figure 11: Band structure calculation for a CoFeB film for a two dimensional square lattice. Top: Free spin- 
wave band structure with first and second band marked in red and orange, respectively. Bottom: Solved by 
a set of Bloch states for a film thickness of t = 50 nm, hole distance a = 3.5 /im and hole diameter d = 1 jam. 
The splitting at the high symmetry points X, X' at the zone boundary is marked with the shaded yellow 
area. 
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Figure 12: Distribution of the spin- wave amplitudes at the point of high symmetry X, X' and M (from [35]). 
From the left to the right, one can see mode profiles of the first band at X', F, X and M. Different colors 
denote positive and negative phase. 
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Figure 13: Band structure calculation for a Ni film for a two dimensional square lattice. Top: Free spin 
wave band structure with first and second band marked in red and orange, respectively. Bottom: Solved by 
a set of Bloch states for a film thickness of t = 50 nm, hole distance a = 3.5 /im and hole diameter d = 1 /im. 
The splitting at the high symmetry points X, X' at the zone boundary is marked with the shaded yellow 
area. 
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4-3. Periodic dipolar potential in structured films: micromagnetic simulations 

Apart from analytical descriptions, numerical simulations provide a convenient tool to 
investigate magnonic crystals and predict their spectral response, please refer also to a cur- 
rent review especially devoted to this topic [36j. For the purpose of numerical simulation, 
it is necessary to implement periodicity. The main physical entity to handle in this context 
is the long-range static demagnetizing field, which is of dipolar origin. The easiest way to 
do this is to construct a simulation volume that consists of many unit cells. One can then 
assume that, at least in the central cell, the surroundings approximate infinite periodicity. 
Of course, such an approach demands high computational power, especially if one has to 
respect fundamental length scales arising from the exchange interaction. The discretization 
length should not exceed a few nanometers if dynamic or static phenomena on this scale, 
like exchange dominated spin waves or domain walls, are considered. The result is that 
simulating many unit cells is only applicable when access to a proper computing facility is 
given. Implementing periodic boundary conditions for a single unit cell circumvents this 
problem by using an accordingly modified demagnetization field. Before the integration of 
the LLG-equation, a demagnetization tensor is computed from the geometry of the sam- 
ple [371 EHj- In a single isolated unit cell, strong divergence of the magnetization occurs 
at the boundaries. By assuming an identical magnetization configuration repeated period- 
ically in space, this divergence is removed, and the demagnetization field of a continuous 
and periodic medium is created. In the strict sense, the imposed exact repetition restricts 
dynamic computations to spin waves with k = 2n7ra~^, where n G and a is the lat- 
tice constant. Only by increasing the number of unit cells within the (periodic) simulation 
boundary, a finer grid in the reciprocal space can be probed. For the computation of the 
static equilibrium magnetization configuration or quasi-static simulations of hysteresis loops, 
the described method provides less ambiguous results. Note that the approach can also be 
applied to the dynamical matrix method, which reformulates the linearized LLG-equation 
as an eigenvalue-problem on a finite grid [39] . Apart from the fact that the dynamic matrix 
method only applies to the linear regime of magnetization dynamics, it is practically limited 
to structures not bigger than 100 nm. 



545 Figure 14 shows the calculation of the internal field (i/int,x = ^ext,x + ^demag,x) for a 
CoFeB sample, where the anisotropy field may be set to zero. If not stated otherwise, the 
external field was canted by an angle of 30° out of the film plane, mirroring the experimental 
situation to be discussed in a later section. The pictures show the equilibrium state of a 
unit cell of a square and an unconventional unit cell of a hexagonal antidot structure. The 

550 colors represent the field component parallel to the external field, normalized to the external 
field amplitude. The cells were computed using the software package Nmag, which can 
apply periodic boundary conditions. Calculations of the square lattice show inhomogeneities 
around the antidots, while in the hexagonal simulation, these regions show higher expansion 
due to the higher fraction of holes per unit area. In both simulations, a lattice periodicity 

555 a = 3.5 /im and an antidot diameter d = 1.36 /xm was used, which lead to filling fractions 
of 11.9% and 13.7% in the case of square and hexagonal lattice, respectively. Note that 
concerning low computational power, discrete element sizes in this case are in the range of 
60 nm. Therefore, exchange interactions cannot be discussed by this simulation. Significance 
has been proven by calculations with bigger and smaller elements showing identical results. 
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Figure 14: Calculation of the total internal field for a square (left) and a hexagonal (right) antidot structure 
with an applied field in the x-direction of /io^ext = 100 mT. The colors show the minima in the extension 
of the internal field around each antidot. Please note the slight asymmetry in the square case to rule out 
unit cell artifacts. 



4^4' Micron and nanometer structuring for static and active magnonic crystals 

Thin film technology provides a large variety of techniques to prepare one- and two- 
dimensional structures. The relevant length scales may range from nano- to millimeters, 
depending on the material under investigation. A technique that is becoming increasingly 
important is lithography using a beam of ions. Starting from a continuous film, the mag- 
netic material is sputtered off by a focused beam of gallium ions (focus 10 nm). Scanning 
the beam across the sample provides the flexibility of e-beam lithography without the need 
for further processing of the specimen. This reduces the efforts and the production of precise 
large scale nanostructured samples is possible. An example of a two-dimensional structure 
is given in Fig. 15 showing high precision, long-range order, and in the inset, an enlargement 
of a well defined constituting antidot. Damage and implantation of Ga ions happens only 
close to the antidot site. With this method different lattice geometries can conveniently 
be achieved and it was used for most of the structures shown in the experimental section. 
Generally, the damping or magnetization of the material is found to be unchanged. Other 
standard techniques, such as optical lithography for structures above 1 jam diameters and 
e-beam lithography combined with lift-off or Ar+ ion milling, are also used as methods for 
structurization. In the bottom row of Fig. [15] we give some examples of magnetic super- 
structures. By ion beam modification, the exchange bias can be modified locally, resulting 
in a separate switching of different areas of a material [40j. The varied switching field allows 
an active modification of the superstructure during the experiments. A modified switching 
field can also be realized by combining soft/ hard materials [41j. An example where dy- 
namic experiments on very small magnonic structures have been realized is shown on the 



very right of Fig. [T5| using self-organized or prestructured porous alumina membranes it is 
possible to enter the range below 100 nm [42]. Ferrites like yttrium iron garnet (YIG) are 
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585 particularly well suited for fundamental research. They have a very low Gilbert damping 
(a = 6 X 10~^) and decay lengths are in the millimeter range. Thus typical structures are 
modified on larger length scales \43\. In YIG already dynamic magnonic crystals have been 
realized by using a set of parallel, periodically spaced, current conducting stripes close to 
the YIG film surface. When the current is turned on the current fiow causes a sine-like 
variation of the film's internal magnetic field |44l [45] , 




Figure 15: Magnetic periodic structures on different length scales and by different methods. Top: Antidot 
structure prepared by focused ion beam etching. The parameters d (diameter of the hole) and a (antidot 
distance) define the filling fraction. Bottom: other techniques involve ion beam bombardment of exchange 
biased layers, which modifies the switching field [40 (left) or combined soft/hard composite materials [Jl] 
(middle), in both cases the magnetic contrast has been imaged. Smallest structures on large areas can 
be produced by interference lithography or by self-organized or prestructured porous alumina membranes 
entering into the nm range [42 (right). 

590 

5. Experiments with magnonic structures 

Magnonic crystals with diflFerent dimensionalities have been experimentally investigated 
in the past. In one dimension, frequency band filters were demonstrated by the group of 
Hillebrands (see review on YIG magnonics [43]). For a two-dimensional magnonic crys- 

595 tal, the simplest case is a squared antidot lattice in a thin ferromagnetic film, as discussed 
by Neusser and Grundler [46j, however the expression 'magnonics' had been formed by 
Kruglyak and Hicken [47j. In the experiment, the control over the internal field distribution 
by the filling fraction (area of the holes per total area) and lattice parameter turns out to 
be crucial. From a broad continuum of spin waves at high energies in low damped CoFeB 

600 films (Gilbert damping a = 0.006), in which the spin- wave propagation length is more than 
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100 /im, Bloch modes appear only for low filling fractions [48]. The static dipolar fields stem- 
ming from a single hole defect have also been studied by combining MFM measurements 
and micromagnetic simulations [49j . One can access the confinement distance and determine 
the interaction of single defects with radius r. It is found empirically that if the separation 
is smaller than lOr, collective eflFects emerge. This is also true for the dynamic case [50]. 
In this experimental section a somewhat larger part will be devoted to the application of 
femtosecond lasers to study magnonic materials. We will divide the presentation of experi- 
mental results into those, where spin-wave localization dominates the collective behavior in 
the periodically structured material and those, where delocalized Bloch modes predominate. 
To give an introduction to this rich field of research, experimental findings from the vari- 
ous magnonic systems obtained with different approaches will be discussed. First a short 
overview on the experimental techniques to detect and excite spin waves will be given. 

5 A. Techniques 

Experimentally spin-wave modes in periodically structured materials have been stud- 
ied by different techniques: microwave based techniques (ferromagnetic resonance (FMR), 
vector-FMR or pulse-inductive microwave magnetometer (PIMM)), in combination with spa- 
tial resolved Kerr microscopy; also, optical pump-probe techniques using femtosecond lasers 
have been applied. Brillouin light scattering (BLS) takes a special role since it also allows 
a resolution in fc-space, and thus, to measure dispersions and band structures directly. The 
different methods have different advantages. Typical specifications for a study of magnonic 
materials are that small areas have to be measured for all in common. Structurization of 
materials in a top-down approach generally allows, as outlined in the previous chapter, only 
small areas. 



5.1.1. Microwave techniques: resonance and time-resolved experiments 

Microwave-based techniques owing to high sensitivity on small areas have the strip 
line geometry in common. Two methods have to be distinguished: ferromagnetic reso- 
nance (FMR), or vector network analyzer-FMR, which use a harmonic resonant excitation 
and methods that are using pulsed activation as pulse-inductive microwave magnetome- 
ter (PIMM). In the first case, a harmonic microwave is guided through the strip line and 
the transmitted signal is measured [51j, similar to a standard FMR experiment using a 
microwave hollow conductor measuring the absorption power of the sample in the cavity, 
however much more sensitive. More complex setups use a local excitation and detection. 
That enables guiding the microwave signal through the magnonic structure directly and 
early studies showed the formation of narrow transmission bands in YIG [52j. Also, this 
allows studying the transmission and refiection of a magnonic material as a function of the 
frequency by a vector network analyzer. PIMM uses a pulsed excitation [53j. The ringing 
down of the magnetic excitations is shown in Fig. [16] for different field values. The signal is 
the inductance induced voltage into the stripline. To fully exploit the advantage in probing 
small volumes, the magnetic structure has to be lithographically structured on top of the 
wave guide. A local approach is possible by using an inductive magneto-dynamic probe that 
can be scanned across the sample, in which spin waves are excited by a microwave strip 
line [54J. 
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5.1.2. Microscopy techniques: imaging of resonant modes 

To map magnetic modes with different frequency spatially, microscopy methods have to 
be used. Experiments on magnonic antidot structures have been performed by Pechan et 
al [55]. They combined Kerr microscopy with time resolution to map modes at different 
frequencies. A current pulse excites the magnetization dynamics, which are probed spatially 
by rastering the sample and measuring the time resolved Kerr signal. The pump beam and 
probe beam come from a femtosecond pulsed laser. The response on different spot positions 



on the sample is reprinted from their work in Fig. [T7| To analyze the data, the time-domain 
data is Fourier transformed to the power spectrum at each spot to a map power distribution, 
which determines the mode location at the resonance frequencies. A very recent work shows 
the power of this method to map the localizing of spin- wave modes [56]. Examples will be 
discussed in section 15.1.51 

5.1.3. Brillouin light scattering (BLS): micro and nano BLS 

BLS can access a further degree of information: the measuring of the inelastically scat- 
tered light under a certain angle allows the calculation of the fc- vector and, by measuring 
the energy shift, the frequency is accessible simultaneously. This method has been devel- 
oped and applied to micron sized structures by the Hillebrands group to highest perfection 
over the last years. Using methods that allow the determination of the resonance of a mag- 
netic mode alone gives access to the energy of the mode for different field directions but 
no information on the wave length. Then, by studying the field dependence, modeling the 
dispersion and comparison to the experimental data, modes can be identified. Naturally 
for that reason, BLS is the most powerful method to access the dispersion and spin-wave 
band structure of a material directly. For example, a wide range of wave vectors can be 
probed [57j , and three-magnon process investigated [58j . The method can be combined with 
microwave excitation to select or populate well defined modes or to investigate the propa- 
gation of spin waves along a magnetic wave guide stripe. Two recent developments should 
be mentioned: by using phase information, phase fronts can be investigated [59]. Also, the 
resolution of the method has been considerably increased just lately: using high resolution 
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Figure 17: Time and spacially resolved MOKE [55]. 

scanning microscopy in combination with BLS, the resolution has been increased, approach- 
ing 200 nm (/x-BLS). This allowed the mapping of modes in permalloy wave guides, shown 
in Fig. [T] A further step has been published by the Demokritov group. Using an optical 
aperture on a cantilever, near field imaging of the BLS intensity allows mapping of the edge 



modes of a micron sized ellipse, shown in Fig. 18, with a resolution of even below 55 nm. 
The /x-BLS technique was also used for the first observation of Bose Einstein condensates of 
magnons j60] under continuous pumping of spin waves at microwave frequency. 



5.1.4' Femtosecond laser techniques: photo-magnonics 

A technique to investigate ultra fast processes involving spin waves is the excitation with 
femtosecond (fs) laser pulses, a field recently reviewed by Kirilyuk et al. [62j. Such ultra 
short pulses can be achieved using titanium doped sapphire seed lasers. Due to the broad 
emission spectrum (FWHM = 200 nm) of the crystal, by mode locking of more than 10^ 
modes fs laser pulses can be generated. Resulting pulses with a repetition rate of typically 
80 MHz have an average energy of about 6nJ. At the expense of repetition frequency, 
amplification of single pulses to ^ 5 /aJ then at 250 kHz is feasible to achieve large degrees 
of demagnetization upon absorption of the pulse by the sample. Focusing the pulses to a 
60 /xm spot (in diameter) onto the sample, one achieves a fiuence of 60mJcm~^. 

However, the use in magnonics is relatively new and we will therefore go into more 
detail. For all-optical experiments, the laser beam is split into two pulses. The stronger 
pump pulse (95% of the intensity) yields demagnetization in the order of up to 50% and is 
used to excite the spin waves in the sample. The other, less intense, pulse locally probes 
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Figure 18: Left: Schematic of how to achieve nanometer resolution with a BLS setup (from [61 ). An AFM 
tip with a nanosize aperture is used for near-field imaging. Right: Experimental data obtained from an 
elliptical permalloy disc showing precessional modes localized at the edges. The applied field tunes the 
spatial extension of the magnetic excitation in a manner proportional to 1/i^ext- 

the magnetization of the sample. Due to the magneto-optical Kerr effect, upon reflection 
from the surface, the polarization changes depending on the magnetization. It is remarked 
that in this technique not the absolute, but only the change of magnetization induced by 

695 the pump pulses is recorded. In the experimental setup [63j the pump pulse is guided over 
a delay stage in order to tune the time delay between both pulses. The magnetic field of 
up to /io^ext = 150 mT is in the longitudinal Kerr-effect configuration, but tilted 30° out of 
the film plane. Due to the very small signal, a double-modulation technique is used with a 
mechanical chopper operating at 800 Hz in the pump beam and a photo elastic modulator 

700 (PEM) at 50 kHz in the probe beam. In the experiments, the time delay between pump and 
probe can be varied in steps of very few fs up to the maximum range of 1 ns after excitation. 
A data set of such a measurement, recorded on a 50 nm nickel film, is plotted in Fig. [19] 
(right). Before the pump pulse arrives, the magnetization is in a state of equilibrium. At 
Ar = 0, the pump pulse hits the sample, leading to the demagnetization of the sample, as 

705 shown by the strong change of the magnetization. As one can already see in the raw data, 
the magnetization shows two different modes of precession. After analysis (subtraction 
of the incoherent background and Fourier transformation) these modes can be identified 
as the uniform Kittel and the perpendicular standing spin-wave modes. Both precession 
amplitude and frequency increase linearly with the external field, as expected. For a better 

710 visualization, all frequency data is also plotted in a color code. 

Compared to BLS measurements, the pump-probe technique only probes the Kerr ro- 
tation, and thus, is not directly /c-sensitive. Comparing it to the resonant excitation at a 
fixed frequency in FMR, significant differences are evident. The excitation mechanism using 
a pump pulse does not select a certain mode. Instead, high energy spin waves with a broad 

715 distribution of wave vectors are excited. During the first few picoseconds after excitation, 
these relax to lower energy states [28j , cool down and populate modes given by the sample 
and magnonic crystal itself. All modes observed are not externally selected (for example, 
by an alternating microwave field), but are selected by the sample itself. Hence the re- 
sults of this technique lead to different information on the magnon density than the ones 
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Figure 19: Time resolved pump-probe setup (left) showing intense pump pulses for excitation and the 
detecting probe pulses with tunable time delay At. Analysis steps for a 50 nm Ni film (right): time-resolved 
trace as a function of the pump-probe delay, Fourier transform and color plot of the Fourier power. 

720 described previously. This technique combines photo and spin-wave excitation, hence the 
term photo-magnonics. A recent work by Hamrle showed that both optical techniques, BLS 
and TRMOKE are comparable in their sensitivity [64j . 

5.1.5. Fundamental mechanisms for spin-wave excitation using lasers 

The full relaxation path of spin waves after laser excitation is a research field of its own. 

725 We will therefore briefly describe the ideas in this separate section. Since there is no di- 
rect coupling between the spin degree of freedom and the light field we want to shed some 
light on the excitation mechanism specifically used in the examples to follow. How can spin 
waves be excited and what is the specific mechanism used in photo-magnonics? Just recently 
two parallel investigations, combined theoretical and experimental proof that it is possible 

730 to quantitatively describe femtosecond demagnetization dynamics in its absolute value and 
characteristic time scales j65l [Ml Elj . A slowing down due to intrinsic fluctuations char- 
acteristic to the exchange coupled spin system can be mirrored by the models found. It is 
found for strongest demagnetization or approaching the Curie temperature. Atxitia et al. 
[67j used a stochastic equation of motion that is mapped onto a macrospin description, the 

735 Landau-Lifshitz-Bloch equation. From spatially resolved micromagnetic simulations, one 
gets a microscopic insight into the origin of the delayed remagnetization. Short wave length 
spin waves and higher order spin-wave modes with a high spatial frequency dominate the 
dynamics within the first ps [28[ l68] . These excitations are distributed very efficiently from 
the hot-spin region into the magnetically not disturbed region. For small excitation spot 

740 sizes, the magnetization can be restored very fast by the emission of spin- wave packets and 
high-energy magnons [69j. High-energy magnons have been studied in detailed experimen- 
tally challenging neutron diffraction experiments in ferromagnetic metals in the 80s [29j, 
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investigating the broadening of the spin- wave dispersion curves at high energies as predicted 
also at that time arising from the strong interaction with the St oner band, caUed Landau 
745 damping: near the surface of the BriUouin zone, it results in a severe energy dissipation and 
extremely short lifetime of the high-energy spin waves. 
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Figure 20: Spin-wave relaxation in all-optical experiments after laser excitation. Micromagnetic simulation 
(left, from [28]). Fourier power spectra for different pumping fluencies 14, 57 and 71mJcm~^, Ni film 150 
nm thickness (right). A Damon-Eshbach mode with fixed wave length A from 0.5 - 6 /im as determined by 
the Ni thickness is found. Precession amplitudes can be up to 5°. A summary of the mode and wave-length 
dependence is given below as a function of the Ni thickness (for a magnetic field of 50 mT) (from [25]). 

In the last years, using spin-resolved electron energy loss spectroscopy [70] and theoreti- 
cal investigation [26] these excitations have been related to spin-wave packets with lifetimes 
below 1 ps and localized within a few nanometers only. The starting point of the relaxation 
chain is an energy transfer from the high-energy spin-wave excitations resulting from local 
spin scattering processes to the lower-order spin-wave excitations. The high-energy excita- 
tions are overpopulated in an ultrafast demagnetization experiment. Furthermore they are 
densely populated (demagnetization typically much larger than 1%). It seems plausible that 
the Stoner excitations, that are populated by the elementary spin-flip relaxation processes of 
the hot electron system (exchange scattering for hot electrons [261 [27lJ and Elliott-Yafet spin- 
flip processes for lower electron energies [TU [721 [73], O]), decay into short wave length spin 
excitations and gradually relax to the lower spatial frequency excitations discussed in [28j . 
A very recent work by Schmidt et al. nicely proves these ideas on the femtosecond time 
scales experimentally [75j. The energy is transferred from the highest excited mode into 
modes owning a lower energy following the spin-wave dispersion. From the corresponding 
Fourier transform, through the magnetization proflles in the micromagnetic model spin- wave 



excitations are accessible (spatial frequency) given in Fig. 20, High spatial frequencies dom- 



inate after excitation at r = 0. During the relaxation process, the center of the spectral 
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weight moves towards lower spatial frequencies. For 6.7 ps the center is at around 0.1 nm~^. 
This corresponds to a spatial spin- wave period of 10 nm. The damping results in an overall 
reduction of the spectral amplitudes, accompanied by a shift of the spectral weight to lower 
frequencies. It has to be mentioned that the dipolar Damon-Eshbach modes are not within 
the simulation window because of the film thickness of only 50 nm. Nevertheless simulations 
are an indication for an energy transfer from high-energy modes to low-energy modes within 
the spin-wave relaxation. Here we see some similarities to the Bose-Einstein condensation 
where spin-waves driven by micro-wave fields relax to their lowest state of energy and over- 
populate this states so that condensation sets in by permanently driving up the chemical 
potential [7^. In these experiments the decay from the GHz pump frequency to the fre- 
quency of the lowest energy level of the condensate can be followed directly. The optically 
induced precession verified by many groups [77| has been successfully used to extract the 
magnetic Gilbert damping for a variety of materials [78j . The amplitudes can be increased 
up to few degrees by increasing the pump fiuence. In the excitation spectra the mode of 
the homogenous precession, standing spin waves of up to fifth order for low damping mate- 
rials [36] and dipolar spin waves with a fixed fc- vector inversely related to the film thickness 
have been investigated [25] . Such a well defined fc- vector resulting from all optical spin- wave 
excitation leads to the idea to study the effect of a periodic modification in order to create 
a magnonic crystal with that periodicity. 

5.2. Localization effects in magnonic crystals 

The idea to excite spin- waves all optically in a magnonic crystal had been put forward by 
Kruglyak et al. [47j. In the following we will study 2D magnonic crystals. And in fact: when 
compared to a continuous film as presented in the previous section, magnetization dynamics 
drastically change after arrays of antidots have been inserted. The two dimensional periodic 
modification of the internal magnetic field has been realized by focused ion beam. Shown in 



Fig. [21] are three data sets, recorded on a continuous 150 nm nickel film, as well as on two 
different antidot lattices of the respective sample. On the structured areas (lattice parameter 
3 fim (3.5 /xm), antidot diameter 0.5 //m (l//m), respectively) the Damon-Eshbach mode as 
previously excited at field values below 100 mT is not observed anymore, instead a new 
magnetic mode is observed. Its frequency does not significantly change with the external 
field. This is a hint that a major role is played by the internal magnetic field which is 
significantly reduced next to the antidots, compared to the applied field. These regions of 
internal fields, however, become smaller with increasing field and finally disappear. The 
interaction between different modes at the antidot sites decreases and they are not coupled 
anymore. This is in accordance with the disappearance of these weakly field dependent 
modes at a certain field value which is found to be different for different antidot diameters. 

From different previous works it is well known, that spin-wave modes in antidot lat- 
tices can have localized or extended character, referring to the spatial behavior of the mode 
amplitude. This was shown in micromagnetic simulations for antidot lattices [79] and in 
earlier experiments by Pechan [55]. The results of the latter using time-resolved Kerr mi- 



croscopy is given in Fig. [22} In 2005, they had already mapped two kinds of modes: the 
spin-wave amplitudes of the high frequency modes have dominantly maxima in between 
the antidots, whereas the low frequency modes have strong intensity around the antidot. 



showing a typical structure that we know already from section |4.3[ it mirrors the areas of a 
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Figure 21: Effect of an antidot lattice structure measured by all-optical spin- wave excitation and detection. 
Left: continuous nickel film of 150 nm. Middle: antidot diameter 0.5 /im. Right: antidot diameter 1 /im. 
With the antidot diameter the mode intensity showing weak field dependence and low frequencies is found 
to increase. 




Figure 22: (a), (b) Modes measured by time-resolved Kerr microscopy (from [55]): the Fourier power at a 
resonance frequency is mapped. High frequency modes are found between the antidot structures (a), low 
frequency modes localized in the low effective field regions around the antidots (b) depict the areas of a 
reduced internal field. In (c), numerical results from ^79^ are reprinted that indicate an analogue situation. 
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Figure 23: Effect of acoustic phonons. (a) Ultrafast heating of the surface of an antidot-structured nickel 
film generates a stress wave which is reflected at the substrate, as seen at a delay of 70 and 140 ps as a sharp 
peak. This film was structured with a = 3.5 /im and d = 1 /im. Because of the low filling fraction the lateral 
heating is homogeneous and no acoustic phonons are excited parallel to the surface which interfere with the 
magneto-optical measurements, (b) In the case of nickel dots on a Si substrate the traced reflectivity changes 
show strong oscillations that can be attributed to a standing acoustic wave governed by the periodicity of 
the structures (reproduced from [82]). 



reduced internal field. Recently, this was investigated in a combined study by vector network 
analyzer FMR, time-resolved Kerr microscopy and BLS [80j. The appearance of delocal- 
ized, extended modes strongly depends on the propagation length of the spin wave excited. 



In the case of Fig. [21} where nickel is the ferromagnet under investigation, the spin-wave 
propagation length is on the order of ten microns due to the rather large intrinsic Gilbert 
damping {a = 0.02). This distance is similar to a few unit cells of the antidot lattice. The 
data in Fig. [21] can thus be explained: arising from the demagnetization field, at the antidot 
edges are potential wells for spin waves with a rapidly changing eflFective field. At these sites 
localized modes are excited while in the regions of constant internal field the uniform Kittel 



mode is present. A schematic of the internal field has already been depicted in Fig.[T4|for the 
overlapping inhomogeneities in periodic lattices. In a time-resolved MOKE experiment, also 
modes extending across several unit cells, hence 'feeling' the periodic potential are induced 



by the structure. These will be presented in detail in section 5.3 (Bloch modes in magnonic 



crystals). The parts (middle) and (right) of Fig. 21 display two major differences. The dif- 
ference of the two antidot arrays being the filling fraction / of 2.2% and 6.4%, respectively, 
the relative amplitude of the localized modes changes when compared to the uniform pre- 
cession. Resulting from the larger perturbance by the lattice in Fig. [21] ( right), the localized 
mode dominates the spectrum over a wider field range than in Fig. [21| (middle). The area 
of inhomogeneous internal field grows with the filling fraction so that the observed increase 
in amplitude further supports the localized interpretation of the non-dispersive modes. 

Another possible explanation for a dominantly field-independent mode is to attribute this 
to a non-magnetic eflFect, the phonon surface wave of the thin metal film. Corresponding 



measurements of the time-resolved refiectivity are shown in Fig. 23 and reveal the propaga- 



tion and interfacial refiection of heat-induced stress waves normal to the film plane described, 
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Figure 24: Effect of the antidot lattice geometry measured by all-optical spin-wave excitation and detection. 
Left: square lattice. Middle: hexagonal lattice. Right: mode spectrum around a single hole structure for 
/io^ext = 30 mT. By scanning the laser spot across a single antidot also here a mode at low frequencies is 
found that solely exists around the antidot. 



which can interfere with the measured transient Kerr rotation [81]. Though small oscilla- 
tions are visible on longer timescales, they do not account for the large amplitudes observed 
and described earlier and are at a different frequency. Additionally, with the velocity of 
sound of the respective materials, the frequency of possibly excited standing phonons can be 
estimated: Assuming the basic mode of phonons with a wave length of Aphonon = 2a = 7 /xm 
and using v^i = 4900 ms~^ as the sound velocity in nickel, one finds an expected phonon 
frequency of 0.7 GHz. Lateral acoustic phonons in the silicon substrate {vsi = 8433 ms~-^) 
have an expected frequency of 1.2 GHz. As a consequence, acoustic phonons as an origin for 
the field-independent modes can be ruled out. The peaks as seen in Fig. 23 (a) arise from 
the stress waves refiected at the film-substrate interface and clearly do not appear in the 
Kerr spectra, which is a sign that there is no crosstalk of the signal here. A reason is that 
the filling fraction is very low, and thus, the whole surface is heated equally. The diflFerent 
possible contributions to the transient signal have been discussed by Miiller et al. [82] and 
reproduced in Fig. [23] (b). In their case, a metallic dot array produced by optical interfer- 
ence lithography was investigated ail-optically; the small dots heat up very fast, while the Si 
substrate keeps cold. The large diflFerence in expansion creates surface acoustic waves with 
large amplitudes that can be matched to the inter-dot distance, which serves as a periodic 
lattice for the surface acoustic waves. This poses a clear diflFerence to antidot experiments. 
It is also possible to examine other, more complex symmetries than the simple square 
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Figure 25: Calculation of the total internal field (left) in a hexagonal structure at /io^ext,x = 90 mT for 
lattice periodicity a = 3.5 /im and antidot diameter d = 1.36 /im. The inset shows the values for a cut along 
the red line. Schematics of the internal field distribution (right) for square lattice with field along (0,1), 
(1,1), tetragonally distorted and for a honeycomb lattice. 



geometry; then the spin-wave modes in the structured media should mirror the lattice sym- 
metries if they are not only determined by the potential around a single antidot, but inter- 
acting with the next neighbor antidot mode. Static measurements of hysteretic properties 
and electronic transport have been performed and resemble the respective lattice symme- 
tries [83j. According dynamic experiments are shown in Fig. 24, They reveal that localized 
modes are also observed on hexagonal lattices (Fig. 24 (middle)). The mode distance in 
frequency of about 1.5 GHz seems very similar for the square and the hexagonal lattice. It 
is solely determined by the antidot diameter. The reader may note the slightly increased 
periodicity in the case of the hexagonal lattice in order to maintain the filling fraction of 
approximately 2%. Also shown Fig. [24] (right) are measurements determined on a single 
antidot. For a fixed external field, the Fourier power is color-coded as a function of position 
and frequency. At the position of the single antidot (marked by an arrow), additional modes 
can be seen. This supports the picture that these low frequency magnetic modes are local- 
ized at single antidot edges. To get a better understanding, we will first discuss the internal 
field distribution around an antidot site in these lattice geometries and then show more 
detailed angular-dependent experiments. In Fig. [25] a calculation of the total internal field 
in a hexagonal structure at /xoi/ext,x = 90 mT for lattice periodicity a = 3.5 /xm is shown. 
The inset reveals that strong changes of the internal field up to distances of 0.5 /xm from the 
dot are found. In the same figure on the right, schematics of the internal field minima for 
different lattices and magnetic field direction are drawn to illustrate a possible interaction 
of modes located in these minima. 

For honeycomb lattices, one expects to find the 6- fold rotational symmetry, which co- 
incides with a 60° repetition of the spin- wave pattern. Here, when compared to hexagonal 
lattices, regions of rather homogeneous internal field are coexistent with closely packed an- 
tidots. Again, when rotating the sample around the film normal, by tilting the structure 
with respect to the applied field, the behavior of the localized modes can be controlled. 



Respective measurements are presented in Fig. [26] and display the expected symmetry. As 
expected the data for the angles 0° and 60°, as well as 30° and 90° show the same features: 
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opposite to the first two 0° and 60°, for 30° and 90° only one dominant mode is found. The 
880 SEM pictures in the insets iUustrate the drop in the internal field around the antidots. By 
rotation of the sample, the overlap between neighboring antidots changes. In the schematic 
inset, it can be seen that as the overlap in the 30° and 90° is larger, the areas of localiza- 
tion come closer. This results in the observation of a higher Fourier power intensity of the 
localized, non-dispersive modes and goes with a reduction of the Kittel amplitude for both 
885 spectra at angles 30° and 90°, respectively. Only one dominating mode is observed all up 
to /xoi/ext,x = 150 mT. Saturation fields of the films are not changed by structuring; they 
still remain at around /io^ext,x = 10 — 15 mT. Simulations for a Py honeycomb lattice have 
been presented . These reveal the effect of the strong interaction between next neighbor 
dots that is responsible for the formation of spin wave channels through the structure. The 
890 dynamic modes map the inhomogeneities of the internal field - which means that they have 
a dominantly localized nature. 
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Figure 26: Effect of the antidot lattice geometry measured by all-optical spin- wave excitation and detection. 
The honeycomb lattice is rotated: for 0° and 60°, the mode splitting in energy is low. A weak localization is 
found. For 30° and 90°, only one mode is found and the mode splitting in energy is high. In the schematic 
picture, a blocking of the Damon-Eshbach wave propagating perpendicular to H is observed (strong overlap 
of the red areas, which represents a drop of the internal field). 

The effects of a changing overlap and crystal symmetry can also be controlled through use 
of elliptical holes instead of circular antidots. The ellipses have an anisotropic dipolar field 
extending into the film for fields applied along the short axis as compared to fields applied 
895 along the long axis. For the modes localized at the field inhomogeneities around the antidots 
this means a different extension to the next unit cells, then 'feeling' the periodic potential and 
symmetry created by the structure. This overlap is responsible for the considerable changes 
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Figure 27: Effect of the form factor in elliptical hole structures in a square lattice. Tuning of the interplay 
between magnetic modes localized at adjacent elliptical antidots can be achieved by making use of the 
structure's anisotropy. Aligning the larger semi-axis with the external field increases the inter-antidot 
coupling which in turn yields a larger localization. The frequency width is found to increase. This effect is 
schematically drawn on the left. Semi-axes are 250 nm and 500 nm, the square lattice has a periodicity of 
3.5 jiim. 



in the relative Fourier power of the localized and uniform modes in Fig. [27| depending on 
the orientation of the external field. 

The structure consists of elliptical holes with semiaxes of 250 nm and 500 nm, respectively, 
arranged on a square lattice with a periodicity of 3.5 /xm. First, one observes a change in 
relative amplitude of the Kittel mode with the angle between external field and lattice. 
Second, two branches of constant frequency occur in Fig. [27| They are separated by several 
GHz and can be shifted by the angle of the applied field. The separation of the modes is 
large for the magnetic field applied along the long axis of the ellipse. Bearing in mind the 
localized character, this behavior can be understood: By tilting the lattice with respect to 
the external field the distance between to localization sites is reduced and the splitting of 
modes can be controlled. This is in close analogy to collective modes as observed in arrays 
of nanomagnets by Kruglyak et al. [85j. The frequency shift controlled by field direction 
and lattice geometry allows to deduce that also the localized modes are more than the single 
mode at one antidot and thus of magnetic origin. Further studies therefore need to include 
a local mode mapping around the antidots to to investigate the nature of these modes 
unambiguously. 

5.3. Block modes in magnonic crystals 

If one wants to study spin-wave materials closer to the "free electron" case two precon- 
ditions must be fulfilled: first, the propagation length of the spin waves must extend across 
multiple unit cells. In terms of a resonance this leads to the formation of a narrow band 
of states in fc-space. Gubbiotti et al. demonstrated in their recent work a magnonic band 
gap spectrum in 2D magnonic crystals with submicrometer periods ^86j. In general, the 
formation of a band depends on the spin- wave propagation. For Ni with a Gilbert damp- 
ing parameter of a = 0.02 the propagation length from phase velocity can be estimated 
to approximately 10 fim. CoFeB has a low damping with a = 0.006. Thus the spin- wave 
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propagation length is larger than 100 /xm. For YIG with a = 6 x 10~^ even millimeters 
can be reached. The second precondition concerns the scattering potential, which must be 
weak. Already a filling fraction of 0.1 means that the holes in the film make 10% of the 
material. However the distortion in the internal field is much larger. It will extend a fac- 
tor two in radius; this variation of the internal field being different from the homogeneous 
case results in an eflFective filling fraction, which is much larger. Very nice examples can 
be found for spin- waves propagating in one dimensional stripes: strong rejection bands are 
formed. The infiuence of a zig-zag configuration of the magnetization results in a period- 
ically alternating magnetization. Topp et al. showed that spin-wave confinements are also 
found [87] . Magnonic gaps in the band structure have been impressively demonstrated in 
a one-dimensional magnonic material by Chumak et al. [88j. Periodic grooves were me- 
chanically drilled into a ferrite-based wave guide. The transmission shows certain frequency 
bands determined by their periodicity, while the groove depth in the material determines 
the transmitted power leading to rejection efficiencies of 30 dB (contrast of 1/1000). Even 
in a standard ferromagnet such as permalloy, by simply increasing the width of the wave 
guide, thus changing the contrast in the periodic energy landscape, a gradual opening of 
a magnonic gap was realized. In the case of 1-dimensional systems the filtering effects in 
width-modulated stripes have been simulated by Kim [89] and subsequently measured by 



Chumak [90j. Both are reprinted in Fig. 28, The distance is in the nanometer range. Clearly, 
this falls into the exchange dominated spin-wave dispersion. It can be nicely seen, that the 
bands calculated in the micromagnetic mode, show a quadratic dispersion, which is very 
similar to a free electron dispersion. Due to the periodic modification of the width, the gaps 
opening up span multiple GHz. The transmission into the structured region is not allowed. 
In Fig. [29] experimental results from Chumak et al. are displayed for a permalloy wave 



guide. Clearly, a drop in transmission not as high as in YIG, but also by a factor of 1/10, 
is observed. 
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Figure 28: Micromagnetic simulation of a one dimensional magnonic crystal: Py stripe with alternating 
width. Forbidden region (magnonic gaps) in frequency result in a frequency filter effect for transmitted 
waves. The transmission into the structured region is not allowed and the intensity in the gap region drops 
down very rapidly (adapted from [89 ). 

The effect of two dimensional structures on propagating modes was studied by Neusser 
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Figure 29: Experimentally measured transmission of a one dimensional magnonic crystal: Py stripe with 
alternating width. A breakdown of the transmitted intensity by more than a factor of 10 is observed at 
8.4 GHz (adapted from j9Q]). 



by vector network analyzer-FMR [79] . We will present a study here as well using all-optical 



spin- wave excitation and detection in a low damped CoFeB film. Figure [30] shows spectra 
of a structured CoFeB film (thickness t = 50 nm), measured by means of TRMOKE. In 
an unstructured film, with the same thickness, one can identify two modes, which are the 
uniform precession mode and the first order standing spin wave. The structure (square 
antidot lattice with lattice constant a = 3.5 //m, field applied along lattice side) alters the 
dispersion, and introduces band gaps at the zone boundary, which is at = ±7r/a. Due 
to a diminished slope of the dispersion at this point, one here finds an increased density 
of states. The fingerprint of this effect can be seen in the measured spectrum in Fig. [30l 
where a new mode appears, which obeys the dispersion of a Damon-Eshbach-mode with 
k = ±7r/a. Since for t <C a the mode profiles for both directions (±90° with respect to the 
field) are essentially constant and the dispersion is degenerate with regard to this directions, 
a superposition of both spin waves yielding a standing wave is most likely. When the field is 



applied in an angle of 45° with respect to the antidot lattice (see Fig. 31), the appropriate 
dispersion for the same propagation direction with respect to lattice yields again k = ±7r/a, 
and is now fourfold degenerate (±45°, ±135° with respect to the field). The degeneracy 
accompanied by an increase of the density of states results in the observation, that only the 
Damon-Eshbach modes are populated after optical excitation. 

These modes are examples for delocalized, extended Bloch-modes, whose excitation and 
detection in non-frequency and non-/c-selective experiments become possible due the artifi- 
cially altered band structure. 

5.4' Magnonic wave guides 

We have discussed the possibility of completely different behavior of photons in a pho- 
tonic wave guide. So called "slow photons" may emerge that propagate in the way of the 
religious procession of Echternach, i.e. three steps forward and two steps backwards, similar 



41 



30 60 90 120 150 30 60 90 120 150 
External field (mT) External field (mT) 



Figure 30: Bloch modes in a CoFeB antidot lattice. In a continuous film (not shown) only the mode of 
uniform precession and a first order standing spin wave (PSSW) can be found. In a structured film an 
additional mode appears (orange line and points), whose wavelength is determined by the wave vector at 
the Brillouin zone boundary k = ir/a. 
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Figure 31: For a magnetic field applied under 45°, Bloch modes in a CoFeB antidot lattice by the wave 
vector at the Brillouin zone boundary k = n/a dominate the spectrum. The antidot lattice is varied as a 
function of the lattice distance for the same filling fraction. Upon reduction of the lattice constant, a mode 
appears that hardly depends on the magnetic field. 



to pictures drawn in the Krauss article [9T]. Effectively, the propagation is slowed down, 
which ultimately is a signature of the periodic potential. Other features like band gap for- 
mation allow for the achievement of photonic resonators with highest Q factors. Moreover, 
wave guides embedded in a structured material allow to facilitate the guidance of photons 
around corners or the splitting of wave packages into two. In the introduction we have al- 
ready discussed the differences between optical and magnonic wave guides: in the latter, the 
internal field may result in a localization of the spin- wave modes. This can be very different 
for magnonic wave guides. Mapping the ideas from photonics to magnonics is not straight- 
forward. Strong differences are found in the anisotropy of the band structure. A Bloch state 
with the period of the lattice will be a standing wave at the zone boundary. However, the 
mode perpendicular to it will not feel the same confinement. By removing one row of the 
holes, the localized state may propagate along this defect, while for the perpendicular state 
the situation remains unchanged. A solution is to apply the field under 45° with respect to 
the defect, as depicted in an artist's view in Fig. 33 (b). Then the states perpendicular and 
parallel to the channel will have the same dispersion and thus show the same standing- wave 



42 




-200 -100 100 200 




Position (pm) Position (pm) Position (pm) Position (pm) 

Figure 32: Magnonic spin-wave guide prepared by FIB: a "line defect". Below: position of the wave guide 
scanning with the pump beam in the reflectivity data. Right: magnetic spectra for different field values 
scanning the position of the probe beam. 



behavior of modes with k at the zone boundary (in close analogy to Fig. 31). Simpler wave 
guides have been discussed in the beginning (Fig. [T]) and even beam-like propagation can 
be realized by spin- wave caustics as reproduced in Fig. 33 (a) [92l|93j. We have performed 
experiments on a "line defect" , which means one missing line in the magnonic crystal. First 



experiments of this kind are shown in Fig. 32, Scanning the probe beam across the defect 



shows magnetic modes inside the magnonic wave guide that are otherwise lower in inten- 
sity in the structured material. While this reveals the feasibility of the experiments, future 
approaches probing the modes propagating in the wave guides need to take into account 
the anisotropic magnon dispersion with respect to the external field. Especially, the mode 
velocity and propagation length as well as the scattering between modes with k parallel or 
antiparallel to Ms and i/gxt will be a focus aim. An absolute advantage of photo-excitation 
and detection of spin waves becomes evident here: excitation and detection can be freely 
moved in the structure by moving the two laser spots (in the experiment 60 fim for the pump 
pulses and 14 /xm for the detection were used here). 



6. Outlook: controlled spin- wave localization 

1005 We have described in our review the possibilities of spin waves in artificially structured 
media: thereby, we started with novel concepts of spin-wave computing; we described the 
peculiarities of spin-wave guides compared to optical wave guides; and then, we switched to 
the effect of structuring materials - what are the conditions and solution of forming novel 
states. To develop grounds for different types of spin- wave materials, we repeated basic 

1010 concepts of solid state physics of band formation of electrons, and clearly distinguished 
between the case of a small perturbation - starting from a freely propagating wave - and 
the tight binding case. We show that similar concepts can be applied to magnonic crystals. 
Experimentally, one finds both localization effects and wave vectors that correspond to a 
Bloch mode at the zone boundary originating from the periodicity of the antidot lattice. 

1015 What kind of development may be projected into the future? One major difference 
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Figure 33: (a) Spin- wave caustics observed with BLS show beam-hke propagation (from [92]). (b) An artist's 
view of a magnonic spin- wave guide: femtosecond laser pulses allow to seed spin waves freely by moving the 
laser spot, hence "photo-magnonics" . Depicted is the situation for the external field applied under an angle 
of 45° with respect to the antidot lattice. 



compared to photonics is that the peculiar shape of the magnetic potential forming at 
each antidot has a rich complexity. While for the photonic case, the index of refraction 
is given by the alternating material itself and can be changed only in small fractions, the 
filling fraction in the magnetic case extends the inter-material border since it is given by 
the dipolar interaction distorting the internal fields around the antidot, which even can be 
changed by rotating the applied magnetic field. The minima in internal fields where spin 
waves localize can be designed by the antidot shape, as we have demonstrated for elliptical 
antidots. By changing the applied field, the potential landscape can be rotated: the overlap 
of the localized states is changed. One can compare this situation with the tight binding 
approach. Each "atom" is defined by the shape of the dipolar fields, which allows certain 
solutions of the spin-wave function. These "atomic" spin functions can then be used as a 
basis set to model the interactions of the emerging localized spin-wave crystal, which may 
have interesting features, to model correlation effects in correlated materials. These ideas 
have been successfully applied to lattices of single atoms in atomic physics. It would be 
interesting to see if such concepts could apply here, as well. For example, the external 
applied field could gradually change the overlap, e.g. "hopping" between the localized spin- 
wave orbital wave functions. However, such ideas are more farsighted. Current experiments 
will have to show that low damping materials will improve the quality of spin- wave crystals 
and formation of Bloch states. Furthermore, different material combinations need to be 
explored to go beyond the antidot lattice case. Spin- wave computation and spin- wave data 
transmission has to be developed to go beyond the proof of principle and to show that indeed 
a high data throughput may be an advantage of spin-wave based computation in the future. 
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